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1 Introduction 



The transport properties of mesoscopic metallic or semiconducting samples 
at low temperature have been shown to exhibit features characteristic of the 
quantum coherence of the electronic wave function along the whole sample, 
see e.g. It is thus interesting to compare the characteristics of mesoscopic 
systems with those for other coherent systems. 

The aim of the present work is to show that if we reduce the dimen- 
sions of a cylinder made of a normal metal or a semiconductor to mesoscopic 
dimensions and if we assume that electrons interact via the magnetostatic 
interaction, a system exhibits coherent properties absent in macroscopic sam- 
ples. 

Coherent systems like superfluids or superconductors can be in general 
described by a two-fluid model namely the total density n is equal to the 
sum of the densities of the coherent ric and normal components 

n = nc + Tin- (1) 

For macroscopic, nonsuperconducting systems (here after refered as "nor- 
mal") TT-c = 0. 

In this paper we discuss the electromagnetic response and its static limit 
of a mesoscopic cylinder of very small thickness d made of a normal metal 
or semiconductor. The circumference and height of the cylinder are denoted 
by La;, Ly respectively {d <^ L^, Ly). 

We argue that such mesoscopic cylinder can be in general described by the 
two-fluid model. The coherent state encountered here, although in principle 
different from the superconducting state, exhibits the following properties 
which bear resemblance to superconductors: 

• free acceleration of part of electrons 

• strong diamagnetic susceptibility at low fields due to the reduction of 
paramagnetic susceptibility 

• dynamic energy gap coming from magnetostatic coupling of electrons 

• persistent currents in the presence of the static flux 
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• flux trapping 



The relation between the induced current density J(g, w) and the electric 
field E{q, w) is generally written as 



J{q,w) 



iK{q, w) 



E{q,w). 



(2) 



w 



It is the form of the kernel K{q, w) that determines the properties of the 



system both static as well as for finite frequencies. The following limits 
imply 



For normal electrons in macroscopic metallic or semiconducting samples both 
limits are zero. 

In this paper we calculate the conductivity a{w) and we study the first 
limit only as we are dealing with a narrow cylinder and we neglect the q 
dependence of the vector potential A. The results for a thick cylinder will 
be presented in a forthcoming paper. 

We show that the relation is fulfilled in the nonsuperconducting meso- 
scopic cylinder and its magnitude depends on the shape of the Fermi Surface 



The reaction of a one-dimensional mesoscopic ring to a static and time 
varying magnetic field and coupled to a thermal bath has been investigated 
by Trivedi et. al. 0. They discussed using the Kubo linear response theory 
the properties of the dynamic and static response function. Such a one chan- 
nel system is a simplification - experimental samples have many transverse 
channels and the question arises how do the system characteristics depend 
on the transverse dimensions. 

We extend the calculations of the conductivity of a one-dimensional dis- 
ordered ring driven by the electromagnetic field to the system of cylinder 
geometry and check whether coherence will still be maintained in a system 
with many transverse channels. 




(3) 
(4) 



(FS). 
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We also discuss the orbital magnetic susceptibility and show that it ex- 
hibits anomalous diamagnetism related to the presence of a dynamic gap. 

We further investigate persistent currents which arise as a result of the 
flux sensitivity of energy levels. We discuss coherence of currents from dif- 
ferent channels and the magnitude of persistent currents for different shapes 
of the FS. 

In one of our recent papers |^ we discussed persistent, paramagnetic 
currents in the free electron model with different shapes of the FS. The FS 
was assumed to be a function of certain parameters which described a change 
of its shape from circular to nearly quadratic. In this paper we use the tight- 
binding model where different shapes of the FS follow in a natural way from 
the crystal structure and band filling [0. We discuss here, contrary to [Q], 
diamagnetic persistent currents. 

Finally we address the issue of whether a system of interacting electrons 
confined to a mesoscopic cylinder can sustain a spontaneous persistent cur- 
rents in the absence of an externally applied flux. These currents correspond 
to flux trapping in the cylinder. 



2 Coherent properties of a mesoscopic cylin- 
der 

Let us consider a system of spinless electrons constrained to move on a hollow 
cylinder with M channels in the height Ly and sites with lattice spacing 
a in the circumference threated by the magnetic flux 0. We assume that 
electrons interact via the long range magnetostatic (curent-curent) interac- 
tion, the interaction is taken here in the mean field approximation (MFA). 
The tight binding Hamiltonian is of the form 

nm 

te " C^^C„+im if^nm+l^nm ^'^nm^nm+l]j (5) 

where t = r,^^ o , comes from the magnetic flux d> 
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e 



eaA ea(j) i t a tR\ 

A is the vector potential in the x direction. The flux (j) is composed of two 
parts 



+ 0/, 



£1(0), 



(7) 



i.e. each electron moves in the external magnetic flux 0e and in the flux 
coming from the current / in the cylinder, C is the selflnductance coefficient. 
The second term in equation comes from the magnetostatic coupling. 
For a clean sample {Vnm = 0) we can diagonalize the Hamiltonian directly 



^xy 2t 



1 — cos k^a 



2n 



hcN 



cos{kya) 



"^kxky 1 



/cx = 0) = — ( s + — 1 , s = 0, ±1, . . . 



ky — ky{r) 



rn 



(M + l)a 



, r = l,...M, 



V N{M + 1) 
I sr) =8sr I sr), 



sin{kyy), 



(8) 

(9) 
(10) 

(11) 
(12) 



2t 



2n 

cos — s - , 



cos 



rir 



M + 1 



(13) 



The eigenstates of equation (^) we denote by | a) , these are the eigenstates of 
the 7i with impurities which can be diagonalized by computer calculations. 
Here we can treat the effects of disorder in rough analogy to the effects of 
temperature assuming that the elastic scattering causes transitions between 
the states of the perfect system We can expand the eigenstates of disor- 
dered cylinder | a) in terms of those of the perfect cylinder | sr), 
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) = E<I^^). al={sT\a), (14) 



a 



n\a) =Ea \ a), po = ——^ TTTT' Po | «) = /a | a), (15) 

where fa = — rwr — vrrr is the FD distribution function, u is the chemical 
potential. 

In this paper we consider the case where the number of electrons A''^ is 
kept fixed and the chemical potential depends on flux; = p,{Nt,4>) can be 
calculated from the equation J2a fa = Nt. 

To introduce relaxation we assume a small coupling exists between the 
cylinder and the thermal bath, with the relaxation rate denoted by 7. Then 
following IQ, by use of the kinetic equation for the density matrix p{t) we 
calculate the response of the system in the relaxation time approximation. 
For details of the formalism see 

Now we apply a small time dependent flux 6<f) = Lx6A, where 6A has 
only the tangential component. Denoting by STi. the corresponding change 
in the Hamiltonian we can write 

en.J-^i.^-'-^i'pXD, (16) 

c c 2 \ c J 

where li = Ip + Id, h is the induced current operator. 

iv = — r E(e"'""^'c++i^c„„ + h.c). (17) 

Ip is the operator of paramagnetic current, running along the cylinder cir- 
cumference. Id is the diamagnetic current operator 



D, D = —- E(e'®""^^4+imCnm + h.c). (18) 



Calculating by the second order perturbation theory the changes in the en- 
ergy caused by STi and comparing them with the Taylor expansion of the 
energy we obtain important relations 
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{a \ Ip \ a) = —c-^, (19) 



Calculating the current induced by the small time dependent flux Jj = Tr/j/3 
we can flnd the conductivity a{w) 



where 



J{w) = ^Iiiw) = a{w)E{w), (21) 



77/? 

E{w) = 6(j). (22) 



Making use of (|T9|) and ( pOD we get after some algebra 



wd(l) - iw ^ d(j) dcj) - £fs £a - ^f3 - u! - i-f 

(23) 

where / is the equilibrium persistent current 

/ = Trpo/p=-cE/.^, (24) 

V denotes the volume of the sample. 

Up to now the formal calculations went in an analogous way as in [0], 
the new physics appears as we perform the summation about the transverse 
channels. The flrst two terms in ( p3|) are present only in mesoscopic multiply 
connected structures where the energy levels are flux sensitive. The flrst 
term, being purely imaginary, determines coherent response of the system. 
The second term involves the intralevel scattering and depends essentially 
on 7, the third term involves interlevel transitions and is connected with the 
elastic scattering. In the absence of the elastic scattering Ip is diagonal in 
the unperturbated basis | sr) and the third term is zero. The ac conductivity 
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in a multichannel metallic ring in the diffusive regime has been extensively 
discussed in 0. 

In this paper we will focus our attention on the reactive response in the 
static limit. We want to establish the conditions under which the system 
exhibits strong coherent behaviour. 

Calculating the w = limit of equation (|23|) we find 

c dl 

lim wlma(w) = — — TTT- (25) 

Thus the conductivity exhibits an imaginary part which in the low frequency 
limit is proportional to the flux derivative of persistent current. 

It follows from the basic theory of conductivity [0 that coherent electrons 
i.e. those which run without scattering obey the following relation 

lim wlma(w) = — — , (26) 
w^o ^ ' rrie 

where Uc is the density of coherent electrons. From equations (0), (H), (pl|) 
and (^) we find 

c dl 

YimYimK{q,w) = = -- — , (27) 

or from (p6| 



X = — . (28) 

rUe 

The finite limit of the kernel K{q, w) corresponds to free acceleration of the 
part (= ric) of electrons - the feature characteristic of superconductors. 
Using the relation between a{w) and x{^) H 

a{w) = -^XH (29) 
we can calculate the static orbital susceptibility 

2 

X = lim x{w) = lim iwa{w) = — lim wlma{w) = — , (30) 

w^O to— >0 w—^0 ui^ 

where relation (^) has been used. 



8 



We would like to stress that the equilibrium susceptibility and the zero 
frequency limit of a dynamic response coincide as we are working with the 
constant number of particles. From equations (^), (0) and (pT]) we obtain 



X = Xd + Xp, 



where 



(31) 



1 

Xd = -Tj^fa{a\D\a) = n, 

Xd is the diamagnetic susceptibility, n is the electron density 



Xp 



V 



^ 2^ /" f _ f ^ ^ 2^ P^M PlA 



Xp is the paramagnetic susceptibility. 

Equations (^) - (|32D give the following relation for Xp 



(32) 



(33) 



Xp 



(34) 



We see that rin is the density of normal electrons which undergo different 
kinds of scattering processes. Both terms in Xp depend on the shape of the 
FS. 

For a macroscopic "normal" sample the distance between energy levels 
is small and the occupation probability is then a smooth function of the 



energy. If is then justified to replace the sum in equation (|2^) by an integral 
what leads to vanishing currents and to zero limit in equation (p7|). Thus 
in macroscopic samples Xp ^-nd xd nearly cancel each other and the system 
displays only a small residual diamagnetism. 

However the mesoscopic cylinder with ~ l/xm has large energy spac- 
ings in x direction and the occupation probability near FS changes signifi- 
cantly over these spacings. The sum can not be then replaced by the integral 
and we get a finite current / which is persistent at low T because of Quantum 
Size Energy Gap (QSEG) 
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The paramagnetic susceptibility is drasti- 
cally reduced due to the presence of the energy gaps and the system displays 
an anomalous diamagnetism. This problem will be further discussed in chap- 
ter 3. We see that in mesoscopic system the paramagnetic and diamagnetic 
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susceptibilities fail to cancel (n 7^ n„) and we can describe the system by the 
two-fluid model. 

The density of coherent electrons Uc is proportional to the flux derivative 
of persistent current. 

We shall discuss now the current / in more detail. It's magnitude de- 
pends on the strength of the phase correlations between currents of different 
channels. 

A large phase correlation among channel currents means that the increase 
of the flux (f) results in an almost simultaneous cross of the FS by the large 
number of channels. The most favorable situation takes place if the separa- 
tion between the last occupied level and the FS from channel to channel is 
nearly the same. There exists then a perfect correlation among the channel 
currents because the M levels cross the FS simultaneously while the flux 
is changed by one fluxoid - we get then the largest amplitude of the total 
current. 

Let us assume that our mesoscopic cylinder contains a small number of 
impurities. The average current, where the average is taken over impurity 
configurations, can be calculated in the linear response approach. For the M 



channel system one obtains |T^ 



/ T M ±00 

/(0)^exp(--^)E E UL. (35) 



where 



eh 2ti I I (j) he 

Is = — sm— s = 0,±l,... 0=—, 00 = —, (36) 

A is the mean free path. For a sufficiently clean material A can be of the order 
of a few microns and the impurities do not decrease the current significantly. 

The current /(</)) is periodic in 0' with period 1 and can be expressed as 
a Fourier sum 



/(0) = exp 



\ i£, ^ AT 2et exp (- 
~2Aj^,^;^Mi_expf 



21T \ 



sm 



F{m), (37) 
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where A is the mean free path, T* is characteristic temperature set by the 
level spacing at the Fermi surface for an electron moving in the x direction, 

F(m) = sm{kp^{m)a) cos{lNkp^{m)a), (38) 

F(m) is the factor depending on a shape of the FS, kp^ is calculated from 
the equation for the FS. 

In the tight-binding approximation the shape of the FS depends on the 
lattice type and on the filling factor [Q. 

For a 2D square lattice for the half-filled band the equation for the FS 
reads 

cos{akp^) + cos{akpy) = (39) 
and the equation for F[m) takes the form 



F{rn) = sin(arccos(— cos(aA;Fj^(m)))) cos(/A^arccos(— cos(aA;i?j^(m)))). (40) 

The FS is than quadratic with the diagonals along k^, ky axes. 

For the filling factor much less than a half we can expand the cosines in 
the dispersion relation (|TB|) for small A;-values and we obtain 



2t 



1 /27r 

2 V iV 



s — 



1 

2 VM 



rir 



a I 



(41) 



The FS is then circular and equation 



reads 



F{m) = sin 



V 



\ 2 







V 

COS 


2/A^. 







^ _ ( akF^{my 



(42) 



We can also imagine a cylinder made of a set quasi ID mesoscopic rings 
stacked along y axis. The FS will be then nearly flat perpendicular to k^ 
direction and the formula for the current is then of the form 
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( 4:T 2et exp(-^) f27Tl(f)\ 
lU) = exp ^ My ^ / \ sin x 

X sm{kF^a) cos{lNkp^a). (43) 

In order to avoid perfect nesting in all figures we present the results for 
slightly less than half filled bands instead of half filled ones. The Fermi 
Surfaces for such cases will have then rounded corners. 

In Fig.l we present persistent currents for four different shapes of 2D FS Insert Fig.l 
(inserted figure). We see that the amplitude of the current increases with 
increasing departure of the FS from circular because the phase correlation of 
the channel currents increases. 

We have found that for a 2D squared FS, drawn by a solid line in the 
inserted figure, maximal interchannel phase correlations exists for = 2 + ^, 
where n is a positive integer. For other values of the ratio the current is 
much weaker. 

The FS of the shapes as in the inserted figure are frequently met in High 



Tc Superconductors with 2D conduction [|^, |TJ . 

For a 3D case with P channels in the cylinder thickness the current (p^) 
takes a form 



( \ ^ ^ ^ 4T 2et exp(-^) . /2^l<p\ 

(44) 

where F{m,p) is a 3D factor depending on the shape of the FS. 

For a 3D cubic lattice for the half-filled band the equation for the FS is 
of the form 

cos(aA;F^) + cos^akpy) + cos^akpj = 1 (45) 
and we get the equation for F{m,p) 

F{m,p) = sin(arccos(l — cos{akpy{m)) — cos {akp^{p)))) x 

X cos(/A^arccos(l — cos^akpy^m)) — cos{akp^{p)))). (46) 
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The FS is then an octahedron with the diagonals along k^, ky and kz axes. 

For the filling factor much less than a half we can expand the cosines in 
relation (HSl) for small fc-values and we obtain 



1 - 



{akp^] 



{aki 



(47) 



and F{m,p) takes a form 



F{m,p) = sin 



V 



X cos 21N 



\ 



'akFJm)Y f akp^ip)"]^ 



(48) 



The FS is then spherical. 

For a body-centered tetragonal lattice and the half-filled band the equa- 
tion for the FS takes a form 



akf 



' aki 



' ckf 



cos 



COS 



COS 



and F{m,p) is 



(49) 



F{m,p) = sin 



2 arccos ■ 



V 



COS 



akp (m) 



COS 



X 



X COS 



arccos ■ 



V 



cosl^^^Uos^^^ 



The FS is then a cuboid. 

For a face-centered lattice we have the FS at half filled band 



(50) 
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and we obtain F{m,p) of the form 




The FS is then a cube. 

Persistent currents can be both paramagnetic and diamagnetic depending 
on the sample dimensions and on the Fermi energy. Paramagnetic currents 
has been studied in details in [Q, in this paper we will discuss mainly the 
diamagnetic solutions. 

In Fig. 2 we present diamagnetic currents for different 3D lattices. Insert Fig.2 

All the above considerations show that the current in a multichannel cylin- 
der depends on the strength of the interchannel correlations and therefore on 
the shape of the FS. 

In the case of the spherical (circular) FS the channel currents add almost 
without phase correlation and the total current is small, whereas for the 
FS being a cuboid (square) with rounded corners the correlation is almost 
perfect and the total current is large. 

It follows from (^Tj) and (|2^) that persistent currents are carried by coher- 
ent electrons. Thus the shape of the FS determines the density of coherent 
electrons ric in the sample and it increases with increasing the curvature of 
the FS. 

3 Self-sustaining currents 

Finally we discuss the possibility of spontaneous self-sustaining currents i.e. 
those which flow without any external field. 

The magnetic flux which drives the persistent current given by equation 
( |37D is the sum of the externally applied flux and the flux from the currents 
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itself. The inductance coefficient L depends on the sample geometry and can 
be large for a cylinder geometry. 

Equations and form two selfconsistent equations for the current 
and it raises the possibility of a spontaneous self-sustaining current at the 
external flux zero. The possibility of the finite spontaneous current depends 
again crucially on the shape of the FS. 

The phenomenon of the self-sustaining currents is a collective effect and 
it requires many electrons to support the current at = 0. Therefore for a 
2D cylinder with c? ~ lA we do not find any spontaneous current solutions. 
To get it we have to consider the cylinder made of a set of 2D concentric 
cylindrical sheets or a 3D cylinder (in both cases we keep d <^ R). 

The graphical solutions of the self-cosistent equations (|^) and ( pTf ) at 0e = 
are presented in figs. 3, 4 for different shapes of the FS. The intersections of 
the two curves marked by circles give the values of self-sustaining currents. 
The temperature at which the transition to the state with such currents 
occurs is denoted by Tc. 

For diamagnetic currents presented in Fig.3 and Fig. 4 the self-sustaining Insert Fig. 3 
currents correspond to flux trapping in the cylinder - the phenomenon known Insert Fig. 4 
in superconductivity. 

It follows from our considerations that the self-sustaining solutions are 
obtained only in samples where the FS has fiat regions. For spherical and 
nearly spherical FS we do not obtain it because the density of coherent 
electrons is too small. 

There are two reasons for the presence of coherent electrons in mesoscopic 
samples. At ffist the QSEG hampers the scattering at low temperatures. At 
second as will be shown below, the selfconsistent flux (pj increases coherence 
in the sample. 

To be mostly transparent let us consider a cylinder made of a set of 
quasi-one dimensional mesoscopic rings stacked along certain axis for the 
filling factor much less than 1. The energy spectrum (cf. (|4l|) ) of a single ring 
is given by the formula 

2h^7v^ , , , CI 

= — ry(^-'^) s = o,±i,... =</>e + — • (53) 

Assuming that each ring posses an odd number of electrons we can calculate 
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the energy gap at the FS, A = Sgp+i — Sgp- We find 



A = Ao f 1 - 20l + 2^^] , (54) 



where Aq = is the QSEG. 

We see that A contains a term A^: 

Ad = Ao^^J-^, (55) 

00 

coming from the magnetostatic interactions among electrons - A^ is the dy- 
namic part of the energy gap, which has to be calculated in a self-consistent 
way. 

If we calculate the energy of electrons in the cylinder made of a set of M 
rings for (/)' < ^ we find 

ot2 2 ±SF 

Ei<p') = M'-^ ^ (.-0T< Wl), (56) 

^e-l^x s=0,±l 

because = 0g + ^ < 0^ as the current is diamagnetic for 0' < |. 



The gain in the energy in equation (|56|) may be called the condensation 
energy due to orbital magnetic interactions. 

In case of an even number of electrons in each ring, spontaneous fiux (pgp 
is created in a system . (pgp shifts the energy levels and as a result a gap 
appears A = S^sp - Ssp- 

One finds 



A = Ad = Ao ■ 2(/)lp, (l)sp = CIsp, (57) 

Isp is the spontaneous current. This gap is dynamic, because it results from 
the collective action of all electrons which produces spontaneous fiux in order 
to minimize the energy. 

In both cases the fiux coming from the currents produces a dynamic gap 
and therefore increases coherence of electrons. This mechanism of gaining 
the energy is valid only in mesoscopic systems as A^ for Aq — 0. Similar 
considerations can be performed for electrons moving on a cylinder with the 
conclusion that the orbital magnetic interactions (taken here in the MFA) 
increase coherence. 
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It can easily be seen from equation (pSD that increase of an energy gap 
results in a further depression of Xp because it decreases both terms in the 
formula for paramagnetic susceptibility thus enhancing the coherent response 
of the system. 

This phenomenon is analogous to the reduction of the paramagnetic sus- 
ceptibility in superconductors and in organic molecules due to pair correlation 



4 Conclusions 



It is already well established both theoretically |jTO| and experimentally [|T6 



that persistent currents which "never decay" can flow in mesoscopic metallic 
or semiconducting rings being a manifestation of quantum coherence. Such 
currents were previously attributed solely to superconductors. 

The question arises what are the similarities and differences between the 
properties of mesoscopic systems and superconductors [|1T], 0. Superconduc- 



tivity is a collective phenomenon and follows from an attractive interaction. 
To compare the two phenomena we considered electrons interacting via the 
orbital magnetic long range interaction and moving on a thin-walled hollow 
cylinder. The interaction, taken here in the MFA, means that each electron, 
besides the external flux feels a magnetic flux 0/ coming from the currents. 
The selfinductance of the one-dimensional loop is negligible however, it 
can be substantial for a system of cylindrical geometry. 

We calculated the frequency dependent conductivity and the electromag- 
netic kernel. The conductivity is strongly related to the presence of persistent 
currents. We have shown that 

dl 

lim lim K{q, w) = -c— (58) 

i/j-»0 q-^O - cl(p 

The finite limit of the kernel is equivalent to the infinite conductivity. It is 
absent in macroscopic normal metals and present in superconductors. Thus 
although both the elastic and inelastic scattering was taken into account, part 
of electrons (= ric) is in a coherent state and moves without dissipation. This 
feature is connected with the multiply connected structure of our cylinder 
where a change of phase produced by the magnetic flux can not be removed 
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by a gauge transformation and results in an equilibrium current / which is 
persistent at low temperatures. 

In a one-dimensional mesoscopic ring the current / is simply a diamag- 
netic or paramagnetic reaction to the external flux 0e and must vanish at 
0e = 0. The situation looks different in a multichannel cylinder. We have 
shown that the self-sustaining, persistent currents can exist at = in rel- 
atively clean samples (ballistic regime) and for the FS which are sufficiently 
flat. The self-sustaining solutions correspond to "orbital ferromagnetism" for 
paramagnetic currents (this case was discussed in 0) and to flux trapping 
for diamagnetic currents - a feature characteristic of superconductivity. 

The coherent response and related quantities depend on the sample di- 
mensions and on the geometry of the FS. We have shown that the magnitude 
of a persistent current and of a dynamic and static response function depend 
on the phase correlation of currents from different channels. This correlation 
increases with increasing the curvature of the FS. The Fermi Surfaces with 
flat regions follow in a natural way in the tight-binding model. 

The coherent response of a mesoscopic hollow cylinder follows from two 
reasons. At first because of finite size the energy spectrum is discrete and 
QSEG causes that the diamagnetic and paramagnetic parts of susceptibihty 
fail to cancel - this is connected with the single electron properties. At second 
the orbital magnetic interaction among electrons leads to appearance of the 
dynamic gap and it further increases coherence. It manifests e.g. in further 
decrease of the paramagnetic susceptibility and can lead to the coherent 
collective phenomena such as "orbital ferromagnetism" or fiux trapping. 

It follows from our considerations that a mesoscopic cylinder can be in 
general described by the two-fiuid model. The coherent behaviour is deter- 
mined here by the interplay between finite size effects and the correlations 
coming from the orbital magnetic interaction. The energy levels are periodic 
functions of the flux with period 0o and the minima of the total energy de- 
termine stable values of flux contained in the cylinder. We will discuss it at 
related problems in more detail in a subsequent paper. 
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Figure captions 

Figure 1. Persistent currents in Iq units (where /q = j^sm^kpa)) vs flux 
for different shapes of the 2D Fermi surfaces for T = 1.21K, Ep = 7eV and 
= 3. In the inserted flgure the shapes of the Fermi surfaces are shown. 

Figure 2. Diamagnetic persistent currents vs flux for different 3D lattices 
and shghtly less than half filled band: cubic (solid line), body-centered 
tetragonal (dotted line), face-centered (dashed line) and cubic for filling fac- 
tor much less than 1 (dash-dotted line) for T = 1.21K,Ep = 7 eV and 
Lj. = 3.14^m, Ly = 0.98/im, = 0.01/xm. 



Figure 3. The graphical solution of a set of self-consistent eqs. (^) and 
(1^ for different temperatures and for body-centered tetragonal lattice for 
Ep = 8 eV and = 3.14//m, Ly = 0.94/im, = 0.02//m. The nonzero 
crossings of the straight line (|^) with the current-fiux characteristic (|37|) de- 
noted by circles correspond to fiux trapped in mesoscopic cylinder. 



Figure 4. The graphical solution of a set of self-consistent eqs. (^) and 
(0) for different temperatures and for face-centered lattice, for Ep = 6 eV 
and = 3.14/im, Ly = 0.96/im, = 0.03/im. The nonzero crossings of the 
straight line (|^) with the current-flux characteristic (37) denoted by circles 
correspond to flux trapped in mesoscopic cylinder. 
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